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Abstract 

The SLIM experiment at the Chacaltaya high altitude laboratory was sensitive to nuclearites and 
Q-balls, which could be present in the cosmic radiation as possible Dark Matter components. It was 
sensitive also to strangelets, i.e. small lumps of Strange Quark Matter predicted at such altitudes 
by various phenomenological models. The analysis of 427 m 2 of Nuclear Track Detectors exposed 
for 4.22 years showed no candidate event. New upper limits on the flux of downgoing nuclearites 
and Q-balls at the 90% C.L. were established. The null result also restricts models for strangelets 
propagation through the Earth atmosphere. 

1 Introduction 

Nuggets of Strange Quark Matter (SQM) composed of approximately the same number of 
up, down and strange quarks could be the true ground state of quantum chromodynamics [1] . 
SQM nuggets could be stable for all baryon numbers in the range between ordinary heavy 
nuclei and neutron stars [A ~ 10 57 ). They may have been produced in the early Universe [T], 
and could be present in the cosmic radiation as a possible component of the galactic cold 
dark matter (with typical velocities of ^lCU 3 c) pj]. As the strange quark is massive compared 
to almost massless up and down quarks, surface tensions lead to the suppression of a few 
s quarks. Thus, SQM should have a relatively small positive electric charge compared to 
that of heavy nuclei [3J |H [5]. In what follows, macroscopic quark nuggets, neutralised by 
captured electrons are called nuclearites. Otherwise, and generally for small baryon numbers 
(A < 10 ), assumed to be quasi totally ionized, they will be called strangelets. In some cases, 
the term SQM or SQM bags will be used interchangeably for nuclearites and strangelets if 
the properties discussed apply to both. In ref. [5! it was shown that an even more stable state 
of nuclear matter could exist, the so called "Colour-Flavour Locked" (CFL) SQM, because 
of the occurrence of a Cooper-like pairing between quarks of different colour and flavour 
quantum numbers. Strangelets are unlikely to have survived from the early Universe but they 
could be produced in the core of neutron stars and released into the Galaxy in very energetic 
astrophysical processes involving strange star collisions @[7] and/or supernovae explosions [8]. 
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Other particles discussed in this paper are Q-balls [S]. They are hypothetical supersymmet- 
ric coherent states of squarks, sleptons and Higgs fields, predicted by minimal supersymmetric 
generalisations of the Standard Model of particle physics. Q-balls could have been copiously 
produced in the early Universe and may have survived till now as a dark matter component. 
Their general properties are different from those of SQM but in some cases, as discussed later, 
the flux limits derived for nuclearites can apply also to them. 

The SLIM (Search for Light Magnetic monopoles) experiment, operated for more than 
four years at the high altitude Chacaltaya laboratory in Bolivia (5230 m a.s.l.), was designed 
to search for light and intermediate mass magnetic monopoles |10) . It was sensitive also to 
SQM and Q-balls [TT] offering the possibility to extend the searches for these particles to a 
lower mass range. 

In the following, we present a short description of the apparatus and of the experimental 
procedure used for the analysis. The properties and energy losses of SQM and Q-balls in 
Nuclear Track Detectors (NTDs) are also discussed. Finally, we give the flux upper limits set 
by SLIM. 

2 Experimental 

The SLIM experiment was a large array (~430 m 2 area) of Nuclear Track Detectors (NTDs) 
deployed at the Chacaltaya laboratory. The detector modules consisted of stacks composed 
of three layers of CR39®, three layers of Makrofol DE®, two layers of Lexan and a 1 mm 
thick aluminium absorber to slow down or stop nuclear recoils. The detectors were exposed 
to the cosmic radiation for an average of 4.22 years [10] after which they were brought back 
to the Bologna Laboratory where they were etched in "strong" and "soft" etching conditions, 
as described in [TTJ [T2] fTS] , and analyzed. 

The strong etching conditions for CR39 were: 8N KOH + 1.5% ethyl alcohol at 75 °C for 
30 hours. The soft etching conditions were : 6N NaOH + 1% ethyl alcohol at 70 °C for 40 
hours. For Makrofol the etching conditions were 6N KOH + 20% ethyl alcohol at 50 °C for 
10 hours. 

We recall that the formation of etch-pit cones ( "tracks" ) in NTDs is regulated by the bulk 
etching rate, vb, and the track etching rate, vt, i.e. the velocities at which the undamaged 
and damaged materials (along the particle trajectory), are etched out. Etch-pit cones are 
formed if vt > vb- The response of the CR39 detectors is given by the etching rate ratio 
p = v T /v B - 

The CR39 and Makrofol nuclear track detectors were calibrated with 158 A GeV In 49+ 
and Pb 82+ beams at the CERN SPS and 1 A GeV Fe 26+ at the Brookhaven National Lab- 
oratory (BNL) Alternating Gradient Synchrotron (AGS). The calibration layout included a 
fragmentation target and CR39 (plus Makrofol) NTD foils in front of and behind the tar- 
get [21 [13] . The detector sheets behind the target detected both primary ions and nuclear 
fragments. After etching, the standard calibration procedure was the following: (i) measure- 
ments of the base area of each beam and fragment tracks in CR39 with an automated image 
analyzer system [TS]; (ii) for each fragment peak the Z//3 was obtained and the reduced etch 
rate (p — 1) was computed. The Restricted Energy Loss (REL) due to ionization and nuclear 
scattering was evaluated, thus obtaining the calibration data expressed as (p — 1) vs REL for 
all detected nuclear fragments [T3] . 

The threshold for CR39 was at Z/f3 ~ 14 (f3 = v/c) in strong etching conditions corre- 
sponding to REL ~ 200 MeV g _1 cm 2 . For soft etching the threshold was at Z//3 ~ 7 which 
corresponds to REL ~ 50 MeV g _1 cm 2 . 

About 50 m 2 of the SLIM CR39 contained 0.1% of DOP additive in order to reduce 
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the neutron induced background by raising the detector threshold [IS]. The threshold of 
CR39(DOP) was at Z//3 ~ 21 in strong etching conditions corresponding to REL ~ 460 MeV 
g" 1 cm 2 , and Z/(3 ~ 13 for soft etching, corresponding to REL ~ 170 MeV g _1 cm 2 . For 
Makrofol the threshold was at Z/f3 ~ 50 (REL - 2.5 GeV g^ 1 cm 2 ). 

The standard scanning procedure for the SLIM experiment was the following. After strong 
etching (applied to all uppermost CR39 layers of each stack) the CR39 foils were scanned twice 
by different operators, using a stereo microscope with a 3x magnification large optical lens. 
The signal looked for was a hole or a biconical track with the two base cone areas equal within 
experimental uncertainties. The double scan guarantees an efficiency of ~100% for finding 
a possible signal. The bottom CR39 sheets were softly etched when a track was considered 
as a possible "candidate" in the uppermost layers (This was about 10% of the cases). The 
third CR39 sheet was etched only in few cases, when there was still a possible uncertainty. A 
few % of Makrofol foils were etched for similar reasons. More details about the experimental 
analysis procedure of NTDs, are given in ref. [13j and references therein. 

No candidate events remained after all cross checks. Two "very unusual events" were 
observed and were finally classified as manufacturing defects in a small subset of CR39 NTDs 
(~ lm 2 ). These "unusual events" are discussed in detail in ref. [17] . 

3 Energy Losses and Detection of Strange Quark Matter 

In this section we shall recall the main properties of nuclearites and strangelets, their energy 
losses in NTDs and some propagation models in the Earth atmosphere. 

3.1 Nuclearites 

Nuclearites are considered to be large strange quark nuggets, with overall neutrality ensured 
by an electron cloud which surrounds the nuclearite core, forming a sort of atom [21 [18] . 
Their radius is never smaller than 10 -8 cm. As the core size increases the electron cloud 
is increasingly contained in the nuclearite core. If the nuclearite mass is larger than 1.5 ng 
(8.4 • 10 14 GeV/c ) all electrons are inside the core. 

In ref. [2] it was postulated that elastic or quasi elastic collisions with atoms and molecules 
of the traversed medium are the relevant energy losses of non-relativistic large mass nuclearites: 



where p is the density of the traversed medium, v the nuclearite velocity and a is its cross 
section: 



Pn = 3.5 • 10 14 g/cm 3 is the nuclearite density and M its mass. 

Fig. QJi shows the REL of nuclearites in CR39 as a function of their initial velocity for 
different masses, computed using Eqs. [T] It can be seen that the REL of any nuclearite 
with velocity (3 > 3 • 10~ 5 is above SLIM's detection threshold. 

A nuclearite with mass M entering the atmosphere with an initial velocity vq << c, after 
crossing a depth L, will be slowed down to 




for M < 1.5 ng 



for M > 1.5 ng 



(2) 
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Figure 1: (a) REL versus j3 for nuclearites of different masses. For comparison the REL of magnetic monopoles 
with charge g = 2go is also shown, (b) REL versus f3 for strangelets with small electric charges and masses 
{Z, A) = (14,300), (17,400), (21,580) and (46, 1000). The thresholds (strong etching) for CR39 and CR39(DOP) 
are indicated by the horizontal lines. 




v(L) = v exp -— / pdx (3) 



where p is the air density at different depths and a the interaction cross section from Eq. [2] 
The accessible region for the SLIM experiment at the Chacaltaya level (540 g/cm 2 of 
residual atmosphere) is expressed as the minimum incident velocity of nuclearites at the top 
of the atmosphere, which yields at the detector level a detectable signal (i.e. REL above the 
detector threshold) as can be seen in Fig. [21 It is assumed, in the propagation through the 
Earth atmosphere, that the variation in mass of the nuclearite and gravitational effects, are 
negligible. For an initial velocity of vq ~ 10 _3 c, SLIM would have detected all nuclearites 
impinging on the detector with a mass equal or greater than 3 • 10 10 GeV/c 2 . 

3.2 Strangelets 

One of the most important properties of SQM bags is their charge to mass ratio, which is 
very small compared to ordinary nuclear matter. For "ordinary" strangelets, as opposed to 
CFL strangelets, the relation charge Z to mass number A is given by [I]: 

Z ~ 0A (iSO^v) A fOT A « 103 

(4) 

Z ~8 f ms V A 1 / 3 for A»10 3 

\150McV/ 



while for Color-Flavor Locked strangelets it is [S] 



4 



mill 




10° 



10' 



10 9 10 



10" 10" 10" 



10 



14 



10 



15 



10 



, 16 



Nuclearite Mass (GeV/c 2 ) 



Figure 2: Minimum velocity for nuclearites at the top of the atmosphere to be detected by the SLIM experiment, 
versus their mass. For (3 m in ~ 10 -3 all nuclearites with mass > 3 • 10 10 GeV/c 2 should be detected. 



where m s is the strange quark mass usually considered to be m s ~150 MeV. The Z/A ratio 
is a key feature for the detection of strangelets. 

Strangelets should be stable for all baryon numbers exceeding a critical value A cr i t ~300 [3]. 
Although, it has been shown that some underlying shell effects may lead to smaller stable 
strangelets [T!5] . They could be accelerated to relativistic velocities by the same astrophysical 
processes responsible for the acceleration of cosmic rays in the Galaxy. If we assume that 
strangelets have no associated electrons then their interaction with matter, in particular in 
NTDs, is expected to be similar to heavy ions but with a different Z/A. The REL of 
strangelets in CR39 nuclear track detectors is illustrated in Fig. [TJd. At low velocities we 
estimated it from "Ziegler's fit" to experimental data from low energy ion measurements |20) . 
At energies > 100 keV/nucleon the REL was computed from fZT\ 



The terms in Eq.[6]are defined in ref. [21]. T upper = MIN(T cnt , T max ) where for CR39 T cut 
= 200 eV and for Lexan and Makrofol T cut = 350 eV. 

From Fig. [TJd, one can see that strangelets with charges above threshold and velocities in 
the range (2 • 10 -5 1) c can be detected with the SLIM CR39 detectors. However, since the 
threshold for CR39 (DOP) is higher, the velocity range is reduced to a maximum (3 ~ 0.7 for 
charges lower than Z ~ 18. 




(5) 




(6) 
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Figure 3: The minimum velocity strangelets should have at the top of the atmosphere to overcome the geomag- 
netic cutoff at Chacaltaya versus their mass. The solid line applies to CFL strangelets, the dashed line to ordinary 
ones. The change in slope (at ~ 3 • fO 4 and fO 5 ) is due to the solar modulation dominance over geomagnetic 
cut-off at large A. 



Being electrically charged, strangelets are affected by the geomagnetic field, fn order to 
pass the geomagnetic cutoff at Chacaltaya (i? cll i //~12.5 GV [22]) strangelets must have 
a minimum velocity (at the top of the atmosphere) which depends on the baryon number, 
Fig. [3] Note that for ordinary strangelets this minimum velocity was computed using the 
more accurate formula for the electric charge (see ref. [1]). The above condition does not 
guarantee strangelet detection by ground based cosmic ray experiments as one has also to 
take into account their propagation through the Earth's atmosphere. 

3.2.1 Strangelet propagation models and detection modes for the detection at 
high altitude laboratories 

Two phenomenological models for the propagation of strangelets in the atmosphere were 
considered in this paper. 

According to the model proposed by Wilk et al. [23] , large relativistic strangelets lose mass 
in each interaction with air nuclei as they travel through the Earth's atmosphere, until they 
reach their critical stability mass and begin to evaporate neutrons. The minimum mass at 
the top of the atmosphere for a strangelet to reach the Chacaltaya Laboratory was estimated 
to be Ao~2200 for normal incidence. The lightest strangelet at the detector level (with 
A ^300) would have a charge of the order of Z ~ 14 (Eqs. H] and [5]). This is close to the 
SLIM threshold. The flux estimated in ref. [21] was of the order of ^2 • 1CP 13 cm~ 2 s -1 
sr _1 for initial strangelets of energies above 10 GeV. The expected number of events at the 
Chacaltaya laboratory would have been ~ 6 strangelets per 100 m 2 per year for masses above 
A crit and energies above 1 GeV. This estimate was made assuming that the flux of strangelets 
is Fs(A = A crit ) — 2.4 • 10~ 5 F to t (F tot is the primary cosmic ray flux at the same energy per 
particle) and has an abundance spectrum, for A > A cr n, oc A~ 7 ' 5 . 
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In the propagation model proposed by Banerjee et al. [251 [2S], small stable or metastable 
strangelets (with A ~ 42, 54, 60, 84, 102, . . .) pick-up neutrons and protons as they traverse 
the Earth's atmosphere, increasing in charge and mass at different rates (proton absorption 
being quickly suppressed as the energy decreases). The strangelet would acquire its final 
charge at ~10 km above sea level. It was estimated that an initial relativistic strangelet with 
charge Z = 2 at the top of the atmosphere could attain a final charge of about Z = 14. A 
larger initial charge (> 2) at the top of the atmosphere could lead to larger final charges. 
The main requirement is that for this charge/mass acquisition process to occur the initial 
velocity of the strangelet should be lower than 0.7c. Otherwise, the excitation energy of the 
strangelet, when colliding with atmospheric nuclei, could exceed the binding energy and lead 
to its fragmentation. Banerjee el al. gave an estimate of the flux of strangelets that could 
reach the level of Sandakphu in India (at ~ 3.6 km a.s.l) for an experimental set-up equivalent 
to SLIM [26]. This flux, 5 -j- 10 events per 100 m 2 per year, could reasonably be extended at 
the level of Chacaltaya. 

Another phenomenological model for strangelet propagation in the atmosphere was pro- 
posed recently |27j but it was not considered in this work. 

In a recent paper, the flux above the rigidity cutoff for strangelets to be detected by on- 
board satellite instruments was estimated, assuming that they are produced in strange star 
collisions and then propagate through the Galaxy as ordinary cosmic ray particles [2S]. The 
expected flux is 

F tot ~ 2 • 10 5 m-V^sr- 1 A -° A67 Z~ 12 

xmaxiRsM^Gc]' 1 - 2 ^, (J) 
where, Rgc is the geomagnetic rigidity cutoff, Rsm is the solar modulation cutoff and, 



A = (PsnV' 2 / 0-5cm- 3 \ / M 

V0.0057 \ n J ^lO-^Moyr- 1 

/1000kpc 3 \ /930MeV 
\ V J \ mod 2 

A is highly dependent on the input parameters (the number of compact stars in binary systems, 
the average amount of mass released in strange star collisions, M, compared to the solar mass 
Mq) which are poorly known. (3s n is the speed of a typical supernova shock wave, n the 
average hydrogen number density in the interstellar medium, V the effective galactic volume 
and mo the nucleon mass. 




4 Flux Upper Limits for SQM 

The analysis of 427 m 2 of Nuclear Track Detectors exposed at the Chacaltaya high altitude 
Lab for 4.22 years yielded no candidate events. This null observation allows us to set upper 
limits on the fluxes of nuclearites and strangelets. For downgoing nuclearites with (3 ~ 10 -3 
we obtain a 90% C.L. flux upper limit of ^1.3 • 10 -15 cm~ 2 s _1 sr _1 , as shown in Fig. [U 
This upper limit is also valid for downgoing strangelets with velocities at the detector level, 
in the (3 range 2 ■ 10~ 4 -j- 10 _1 as indicated in Fig. \E\ These velocities are dependent on the 
propagation model and on the atmospheric depth crossed by the initial strangelet. 
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Figure 4: 90% C.L. flux upper limits versus (5 for a downgoing flux of nuclearites obtained by the SLIM experiment 
at Chacaltaya. The solid line is valid for nuclearites with Mn < 8.4- 10 14 GeV/c 2 , the dashed line is for Mn ~ 10 17 
GeV/c 2 . 




Figure 5: 90% C.L. upper limits versus (5 for a flux of light downgoing strangelets obtained by SLIM at Chacaltaya. 
The dotted line is for a charge Z = 14 at the detector, the solid line is for Z — 21 and the dashed line for Z = 46, 
respectively. 
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5 Q-balls 



Q-balls are nontopological solitons predicted by minimal supersymmetric generalizations of 
the Standard Model of particle physics. They are hypothesized supersymmetric coherent 
states of squarks, sleptons and Higgs fields, and may carry some conserved global baryonic 
charge Q. 

According to their properties of interaction with matter [29, 30J , Q-balls are classified into 
two groups : (i) neutral Q-balls, generally called SENS (Supersymmetric Electrically Neutral 
Solitons), that should be massive and may catalyse proton decay; and, (ii) charged Q-balls 
called SECS (Supersymmetric Electrically Charged Solitons), that might be formed by the 
SENS gaining an integer electric charge from proton or nuclei absorption. Coulomb barrier 
may prevent further absorption of nuclei. 

In refs. [351 OS] it was pointed out that neutral Q-balls (SENS) would interact with matter 
mainly by converting nuclei into pions with a cross section determined from their size. This 
process occurring at a higher rate than ordinary collisions of Q-balls with nuclei, may not be 
relevant for detection in NTDs. However, it could lead to SENS acquiring an electric charge, 
thus transforming into SECS. SECS could capture electrons and the resulting system would 
be similar to an atom with an enormous heavy nucleus [18] . The cross section for interaction 
with matter will be determined by the size of the electronic cloud. Thus, charged Q-balls 
should interact with matter in a way not too different from nuclearites, i.e. their energy losses 
would be about the same as for nuclearites with a constant radius ~10 -8 cm (given by the size 
of the electronic cloud). The conclusions obtained for nuclearites can be extended to SECS 
of medium-high masses. For SECS with very large masses, Mq > 10 26 GeV/c 2 or Mq > 10 30 
GeV/c 2 (depending on the energy scale of the SUSY breaking theory) the core radius will be 
larger than 10 -8 cm. Consequently, the electrons are inside the Q-ball and the whole system 
would behave as SENS [32] . 

In the case where SECS would not capture electrons, a different method is needed to 
estimate the rate of energy losses [35]. For typical galactic velocities v ~ 10~ 3 c, there are two 
contributions to their energy loss which come from the interaction with electrons and nuclei 
in matter pS ESI [30] ■ Electronic energy losses are dominant at the highest energies and 
can be computed from a Ziegler fit to the experimental data [20] . Energy losses from nuclear 
collisions are computed from the diffusion cross sections considering a screened electromagnetic 
potential as described in [30] and references therein. The REL for Q-balls of two different 
charges Zq = le and 13e is illustrated in Fig. The CR39 and CR39(DOP) thresholds are 
also shown. It can be seen that these SECS, for velocities /5 > 3.5 • 10 -5 , could be detected 
in SLIM CR39. However, in the case of CR39(DOP) and for SECS with charge Z Q = le, the 
range of detectability is reduced to /? > 4 • 10~ 3 . 

Since no candidate event was found in the SLIM experiment, we set the flux upper limits 
as shown in Fig. [7] for Q-balls with Zq > le. It can be seen that, for Q-balls with Zq ~ 13e 
and (3 > 10~ 4 , the limit is about 1.3 • 10~ 15 cm~ 2 s sr -1 . While, for Q-balls with very low 
charges it is considerably higher. 

6 Conclusions 

The analysis of the CR39 NTDs of the SLIM experiment has resulted in new upper limits (90% 
C.L.) on the flux of SQM nuggets and charged Q-balls at the level of 1.3 • 10 -15 cm~ 2 s~ 1 sr _1 . 

In Fig.|5]the flux upper limits for nuclearites with j3 > 10 -3 set by other searches based on 
NTDs [33[[31[35], together with the SLIM limit and the combined MACRO + SLIM limits 
are shown. The galactic dark matter upper bound is also indicated. Above 4 • 10 22 GeV/c 2 
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Figure 6: Restricted Energy Losses of charged Q-balls with electric charge Zq = le and 13e, as a function of 
their velocity in CR39 nuclear track detectors. The detector threshold is also shown (dotted horizontal lines). 
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Figure 7: 90% C.L. upper limits obtained for a downgoing flux of charged Q-balls (SECS) with Zq = le and 
Zq — 13e, plotted vs /3. 
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Figure 8: Flux upper limits vs mass, for intermediate and high mass nuclearites (with (3 ~ 10 3 ), given 
various searches with NTDs, see text. The combined flux from MACRO and SLIM is also shown. 



the upper limit for a 4-7T acceptance is also drawn. 

Fig. [U shows the flux upper limits for charged Q-balls (for Zq > 10 e and /3~2 • 10 -3 ) 
from various experiments employing NTDs ( 33J and [3H [35] [35] as re-estimated in ref. [32 ) 
together with the SLIM result. 

Finally, the flux upper limits set by SLIM for relativistic strangelets are compared with 
the results from previous searches with experiments onboard stratospheric balloons and in 
space (ARIEL-6 [37], HEAO-3 [35], SkyLab [35] and TREK gD]) as shown in Fig. [TO] The 
candidate "exotic" events reported so far are also shown [4"Tl H"2"I |4"5] . 

In order to compare searches performed by SLIM with other experiments we had to assume 
a particular propagation model in atmosphere. 

According to Wilk [24] , the number of events predicted to reach Chacaltaya and be detected 
with SLIM is ^6 events/100 m 2 /year/sr. This would correspond to ~10 events with Z > 14 
and velocities < 0.7c at the detector level. For higher velocities, the flux would be ~1 
event/100 m 2 /year/sr for strangelets with masses A > 2800 at the top of the atmosphere. 
Taking into account SLIM angular acceptance [T3] , the number of events that could have been 
recorded by SLIM NTDs for strange-let charge Z^18 at the detector level is ~ 5 events. 

Assuming the Banerjee propagation model J25[ [26], the number of events with charge 
Z > 14 that should have been detected at Chacaltaya is ~20 -j- 40 events/year. As our 
detector acceptance is larger at the expected velocities (< 0.7c) and the penetration through 
the atmosphere has no suppression effect in the Banerjee model, SLIM should have recorded 
about 80 events. The null observation tends to exclude the Banerjee propagation model of 
strangelets in the atmosphere, while the Wilk model is only marginally excluded. 

At the moment we do not discuss the predictions made in ref. [28j since the comparison 
may not be straightforward. 

In conclusion, SLIM has extended the search for intermediate mass exotic particles with 
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Figure 9: Flux upper limits of charged Q-balls with Zq > lOe versus mass, see text. 
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Figure 10: Flux upper limits vs A for relativistic strangelets as determined by experiments onboard balloon and 
space and by SLIM at mountain altitude (see text). The lines are the expected fluxes from different models: 
Madsen [28] solid (CFL Strangelets) and dashed (ordinary strangelets) lines; Wilk [24] dotted line. The SLIM 
upper limit is shown for three different velocities (f5) at the detector level. 
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observations at high altitude and lowered, in combination with MACRO, the flux upper limit 
for larger mass exotics. New results on the strangelet flux at lower masses (A < 1000) may 
hopefully come from next generation satellite experiments e.g. AMS-02 [U], PAMELA [4"5] . 
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